> 



On 



QCD vacuum tensor susceptibility and properties of transversely polarized mesons 

A. P. Bakulev and S. V. Mikhailov 
Bogoliubov Theoretical Laboratory, Joint Institute for Nuclear Research, 141980, Dubna, Russia 

(February 1, 2008) 

We re-estimate the tensor susceptibility of QCD vacuum, x> an d to this end, we re-estimate 
the leptonic decay constants for transversely polarized p- , p'- and bi-mesons. The origin of the 
susceptibility is analyzed using duality between p- and 6i- channels in a 2-point correlator of tensor 
currents. We confirm the results of Jlj for the 2-point correlator of tensor currents and disagree with 
Q on both OPE expansion and the value of QCD vacuum tensor susceptibility. Using our value for 
the latter we determine new estimations of nucleon tensor charges related to the first moment of the 



transverse structure function h\ of a nucleon. 
PACS number(s): 11.15.Tk, 12.38. Lg, 14.40.Cs 
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INTRODUCTION 

In this paper, we investigate the low-energy properties of the lightest transversely polarized mesons with quantum 
numbers J PC = l~~(p), 1 H (&i) in the framework of QCD sum rules (SRs) with nonlocal condensates (NLCs) as 

00 ■ 

well as with the standard ones. This work was started in 13] where the "mixed parity" NLC SR for the light-cone 

00 ' . . . 

distribution (LCD) amplitudes of both p- and 6i-mesons was constructed. It was concluded that to obtain a reliable 



result we should reduce model uncertainties due to the nonlocal gluon contribution into the SR for LCD. Different SRs 
for each P-parity could be preferable for this purpose. As a first step, to obtain twist 2 meson LCD, we concentrate 
Oh! on the meson static properties using the "pure parity" NLC SR for each meson separately: 

^H. 1) we re-estimate the leptonic decay constants for transversely polarized p(770), p'(1465)-mesons and the 

^ ■ 6i(1235)-meson (l+~) g; 

2) we correct the previous consideration by Belyaev and Oganesyan (B&O) || and provide a new estimation for 

' the vacuum tensor susceptibility (VTS) introduced in f|||. 
'_ 1 

The static characteristics, the decay constants f 7 ^ and "continuum thresholds" s m (parameters of phenomenological 
models for spectral densities) of the lightest transversely polarized mesons in the channels with J PC — 1 and 1^ , 
are tightly connected with the value of VTS. Namely, the difference of the meson properties in these channels fixes 
the non-zero value of VTS: in a hypothetical Nature, e.g. , where the properties of these mesons are the same, VTS is 
identically equal to zero. For the reason that these meson constants should appear in VTS in a form of a difference, 
we have to define them more precisely and in the framework of a unified approach. 

The method of NLC SRs was successfully applied for the determination of meson dynamic characteristics (LCD 
amplitudes, form factors, see, e.g. , |3]j7j and refs therein). One of its basic components is a non-zero characteristic 
scale, X q , of quark-gluon correlations in the QCD vacuum JtJ. This parameter fixes the average virtuality of vacuum 
quarks which flow through vacuum with the momentum k q , (fc 2 ) — A 2 s» 0.4 — 0.5 GeV 2 jl(|, this value is of an order 
of hadronic scale, m 2 p « 0.6 GeV 2 , and is of importance in the calculations. Note here that quite recently the nonlocal 
character of the quark condensate has been confirmed in the lattice calculations |llf| , where an attempt to estimate 
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X q was made. NLC approach can improve the stability and accuracy of SRs even for the case of decay constant 
determination where the NLC effect is of an order of the radiative correction contribution. Therefore, we revise their 
values in pure parity NLC SRs, despite the presence of different estimations for these quantities in literature |8|,||,[|, 
obtained by different ways. For comparison we also calculate all these quantities by the standard way, that corresponds 
to processing our NLC SR on the limit A 2 — > 0. 

The source of the above-mentioned difference of meson properties is the peculiarity of four-quark condensate contri- 
bution to the "theoretical part" of SRs. This contribution is invariant under the duality]^] transformation in contrast 
to all other condensate contributions which change the sign under the same transformation. This peculiarity of 
four-quark condensate contribution will be considered in detail. 

The plan of presentation is the following: firstly, we discuss the QCD SR approach to investigation of the 4-rank 
tensor 2-point correlator for transversely polarized p-, p'- and &i-mesons. Then, we define the duality transformation 
and draw its consequences for the constructed SRs. Finally, we derive a new estimation for the QCD VTS and nucleon 
tensor charges and discuss what is wrong in the consideration of Q . 

DECAY CONSTANTS OF TRANSVERSELY POLARIZED (J PC = 1", l+-)-MESONS 

We start with a 2-point correlator of tensor currents J^ v (x) = u(x)<j^ v 'd(x) , 

n/wi«0( g ) = f d 4 x e lq - x (0\T[J^ + {x)J af) (0)}\0). (1) 



(Note here that due to isospin symmetry this is the same correlator which was studied in §.) This correlator can be 
decomposed in invariant form factors II±, 

W'^iq) = n_( (/ 2 )Pf l/;a ' 3 + n + (q 2 )P^ a(3 (2) 

where the projectors P\.i arc defined by the expressions 

^ ^ [g^q - g" a q»q P ~ 9^q V q a + g V \»q a \ 5 (3) 

pr aP = \ [g» a g vP - g^g va \ - , (4) 

which obey the projector-type relations 

(Pi ■ Pjf^ = puppet? = 5 ..pW<*P (nQ sum Qver pfiiw = 3 (5) 

Then, for the form factors H±(q 2 ) it is possible to use dispersion representations of the form 

n±(<7 2 ) = - / p± ^ ds + subtractions , (6) 
t Jo s- q z 

which after the Borel transformation (with Borel parameter M 2 ) become 



1 We are grateful to O. V. Teryaev who involved us in an investigation of this transformation and suggested this name for it. 
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1 f 

n ± ( q 2 ) - sn ± (M 2 ) = J P± ( s )e- S / M2 d s . (7) 

A phenomenological model for the spectral density /9 phen (s) is usually taken in the form of "lowest resonances + 
continuum" 

P t en (s) = ±2tt \f?\ 2 s.6(s- ml) + pl ert (8)e(8 - s ± ) , (8) 

where and m m are the decay constants and masses of the lowest meson resonances, m — p, p' , b±, contributing to 
the correlator of interest, and p^ rt (s) are the corresponding spectral densities of perturbative contributions to the 
correlators H±(q 2 )- The decay constants are defined via the parameterization of the unit helicity (|A| = 1) states of 
p- 1 p' - and &i-mesons 

(0\u(x)a flu d(x)\p+(p,X)(p'+)) = ifl pl {e M (p,\)p v -e v (p,\) Pll ) ; (9) 
(0 \u(x)a^d(x)\ b+(p, A)) = / b % ra/ 3£ a (p, \)pp , (10) 

here e^(p, A) is the polarization vector of a meson with momentum p and helicity A. To construct SRs, one should 
calculate OPE of the correlators H±(M 2 ) 

^ ± (M 2 ) = ^ £fi? (s)e-^ 2 ds + ^G 2 ) + ^{^qqf . (11) 

We perform these calculations in the approach of QCD SRs with NLCs (see where the coefficients a±, b± become 
functions a±(M 2 ), b±(M 2 ) of the Borel parameter M 2 which tend to their standard values for large M 2 , M 2 ^> A 2 , 
e -<7- > b± = lim b±(M 2 ). The functions a±(M 2 ) 1 b±(M 2 ) can be considered as accumulating an infinite subset 

of the standard condensate (A 2 /M 2 ) J -contributions in OPE. All needed NLC expressions are given in Appendix 
A, while the standard coefficients a±,b±, corresponding to the limit X 2 /M 2 — > 0, are explicitly written below. Their 
values are in full agreement with the preceding calculations performed in Q| and (s| 

, a s (p 2 ) (7 2, s 



9 ' 3 l0g ^ 



(12) 



jkj a± = k ; (13) 

1 , -16 + 80+ 144 208 . 

- — -6_ = = (14) 

(-2) 81 81 v ' 

1 , -16 + 80- 144 -80 

b + = = ~~oT~ ■ ( 15 ) 



(+2) + 81 81 

Here p is the renormalization scale (p 2 ~ 1 GeV 2 ) and the coefficients listed in the central parts of the last two 
lines correspond to the vector (q^ m q), quark-gluon-quark (qG^q) and the four-quark (qqqq) vacuum condensate 
contributions (see details in Appendix A, Q). We write down these coefficients explicitly in order to reveal the 
discrepancy between our results and those obtained by B&O Q, who found, instead, in the last line 

-16-48- 144 -208 
81 ~~ 81 ' 

a result larger than ours by a factor of 2.6. We conclude that in [0 there is a wrong contribution due to the 
quark-gluon-quark vacuum condensate. 
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Collecting all parts 



f p | m p e 



-ml/M 2 



\fbi\ 



|), together, one obtains the following SRs: 

!/'-// ) = - / ' Po CTt {s)e~ s / M2 ds - ^-( — G 2 
n Jo 

-ml /M _ 



77 Jo 
1 /"* 

I" .A) 



b_(Af 2 
2M 2 



1 pert, s -s/M 2 , , a + i a B n 2\ , M M ) / , \5 

(s)e 7 rfs + -y ( — G ) + 2M 2 ^(V^g) 



(16) 
(17) 



The role of NLC, concentrated in a±,b±(M 2 ), is important here, i.e. , at M 2 = 0.6 GeV 2 the total condensate 
contribution in the SR reduces twice in comparison with the standard (local) one. In accordance with QCD SR 
practice the processing of these NLC SRs are performed within the validity window < M 2 < M? (see details 
in ^,^]). These windows are determined by two conditions: the lower bound Af 2 . by demanding that the relative 
value of (GG)- and (gq)-contributions to OPE series should not be larger than 30%, the upper one by requiring 
that a relative contribution of higher states in the phenomenological part of SR should not be larger than 30%. The 
processing with the standard values of vacuum condensates (see Appendix A) gives the decay constants 



/J = 0.157 ±0.005 GeV, fj, = 0.140 ± 0.005 GeV, s] 



2.8 GeV^ 



f£ = 0.184 ±0.005 GeV, s. 



bi 



(18) 
(19) 



which are presented at normalization point ji 2 = 1 GeV 2 
obtained for all of these quantities. 



2.87 GeV'' , 

Very stable curves in wide validity windows have been 
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FIG. 1. The curves of fj in M 2 ; the solid line corresponds to NLC SR with the p'-meson taken into account, 
the long arrows show its validity window; the short-dashed line corresponds to the standard SR without p'-meson, 
the small arrows show reduced validity window for this case; the dashed line corresponds to B&B analyses. 

The processing of the "local" version (at A 2 — > 0) of the SRs (|l6|)-([l7|) leads to the values 0: 




2 To provide a clear comparison with the results of B&O, who do not take into account p'-meson contribution, condensate 
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fj = 0.179(0.170) ±0.007 GeV, fj, - GeV, sj p , = 2.1 GeV 2 ; (20) 
/£ = 0.191(0.178) ± 0.009 GeV, s£ = 3.2 GeV 2 , (21) 

which accuracy looks worse. Really, the curve corresponding to fj in M 2 is shown in Fig. 1 (solid line) in comparison 
with the result of the standard approach without //-meson (short-dashed line). For the first case, the validity window 
expands in all the region (0.55 — 1.20) GeV 2 (long arrows) while for the latter case it shrinks twice to the region 
denoted in figure by the small arrows M 2 and Ml. Note, that the standard SR "pushes out" the p' meson and 
does not allow to obtain his parameters, while the NLC SR is sensitive to this meson and even allows to determine 
its mass ||. We demonstrate on the same figure the curve for /J (dashed line) obtained in Q by Ball and Braun 
(B&B) in the framework of the standard approach, the same small arrows denoting its real "working window" . Note 
here that processing B&B SR just in this thin working window results in a curve very similar in shape to the upper 
short-dashed one with the average value /J(lGeV 2 ) = 0.171 GeV. Note, that recently performed lattice estimates |l^] 
give /J iQtt (4GeV 2 ) = 0.165(11) GeV, that approximately agree with both the NLC © and the "standard" (§(j) 
values. So, we can conclude that our improved SRs (|l6|)-(|l7|) are really justified and produce reliable and stable 
results. All the results obtained by processing "pure parity" (|ll| [l?]) and "mixed parity" NLC SR || are collected in 
Table 1, in comparison with the previous results in SB]. 





"pure parity" SR 
based on H T (— for p, + for bi) 


"mixed parity" SR 
based on (II_ — H + )/q 2 


Source 


Here 


B&B § 


B&O | 


HereQ 


B&B § 


ff [MeV] 


157(5) 


160(10) 




166(6) 


163(5) 


f pl [MeV] 


140(5) 










s P [GeV 2 ] 


2.8 


1.5 




1.5 


2.1 


fl [MeV] 


184(5) 


180(10) 


178(10) 


179(7) 


180 ««d 


s bl [GeV 2 ] 


2.87 


2.7 


3.0 


2.93 


2.1 



Table 1. Estimates for decay constants / T (lGeV 2 ) of transversely polarized p(770), //(1465)- 
and &i(1235)-mesons based on processing QCD SRs in different approaches. 



nonlocality and a s -corrections in the perturbative spectral density, we write down the results of processing our SRs in the same 
approximation in parentheses. 

3 We are indepted to D. Becirevic, who informed us about these interesting papers, contaning lattice estimates mass and decay 
constants of mesons. 

4 The estimates presented in this column have been obtained by processing the "mixed parity" SR established in |J. We 
improve the model for phenomenological spectral density using the features of phenomenological spectral densities of "pure 
parity" SRs. 
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It is interesting to note that in spite of the discrepancy in the OPE coefficients, the authors of [|| obtain for 
a value of 178 ± 10 MeV which is quite close to the value found by B&B ||: 180 ± 10 MeV. This compensation 
effect happens due to the fact that both groups of authors used different sets of the condensate input-parameters 
in the SR and this resulted in approximately the same overall contributions of the quark condensate: B&B had 
{\b+)Tta s {qq) 2 = -4.22 10~ 4 GeV 6 ; and B&O, (\b + )Tta s {qq) 2 = -4.92 10- 4 GeV 6 , see Appendix B. 

DUALITY AND ITS BREAKDOWN 

Let us consider now an operator D transforming any rank-4 tensor T^ v '' a ^ to another rank-4 tensor yg"> Q = 
(i)Ty v -' a P with 

Vft&p = ^\, v ,e a ,^ and & = 1. (22) 
Our projectors p^ v ' a P and p£ v ' a P under the action of this operator transform into each other 

(dpA = Pp a ? ; (DP 2 j = Pp' a/3 , (23) 
whereas the correlator Il^' a "(q) transforms into the correlator of dual tensor currents J^ u (x) = u(x)a ,iU ^r- ) d(x) 

(Dtt)^ a0 (q) = J d 4 x e lq - x (0\T[J^ + (x)jf(0)}\0) . (24) 
There is a good question: how are H^' a0 (q) and (DH)^' 013 (q) connected? 

In perturbative QCD with massless fermions, taking into account the standard anticommutations, one easily arrives 
at 

(DUr p ^(q) = -U^(q), (25) 

from which it follows that 11^"^ (g) is anti-self-dual. The same (anti-dual) character is inherent in the phenomeno- 
logical models, see Eq. (||). 




FIG. 2. Diagram with insertion of four-quark condensate. 

The same reasoning is valid almost for all OPE diagrams, those with a gluon condensate, with a vector quark 
condensate, and with a quark-gluon-quark condensate. Only the diagram with four-quark scalar condensates is 
different (see Fig. 2): in that case there are 2 7-matrices on one line between two external vertices (1 from the fermion 
propagator and 1 from the quark-gluon vertex) because the scalar condensate cancels one 7-matrix. Thus, we realize 
that the OPE contribution involves two parts, one being anti-self-dual (ASD) and the other one self-dual (SD) 



G 



Tlgpf (g) = ASD A " /;a!/3 (g) + SD"" ;a ' i (g) , (26) 
(-£>ASD)^ ;Q/3 ( g ) = ASD^ ;a/3 (g) = n ASD (g 2 ) f P w°>P _ P g»;°<P\ f (27) 
(£>SD)^ !a ^(g) = SD^ ;aP (q) = n SD (g 2 ) fp^ " 3 + . (28) 

The appearance of the SD-diagrams breaks the anti-duality of the two correlators II and yDTlj . 

We can rewrite formulae (^)-(^) to obtain the following representation for the OPE-induced part of correlator: 



n^(g) = pr aP |n 0i V) + u ASU (q z ) l + pp ap |n 0i V) - n ASi V) l . (29) 

As a simple consequence of this representation and Eq.(j^) we have a useful relation 

n^( g ).6n SD ( ? 2 ) . (30) 



T SD/2s , tjASD/ 2n 



T SD/2s ttASD/„2n 



Using relations (]2j])-(|28|), one can easily calculate the OPE coefficients for different diagrams. For example, let us 
consider the (gGg)-condensate and its contribution to the coefficient b±. Indeed, we know that this contribution is of 
the ASD-type, that is 

n ^> =c{q 2 )(pr aP -pr a0 

Therefore, for a light-like vector z, one has 



TiV;af3 , 2\ 2(<7 ' z ) 2 

(qGq) = ^G^aiMxZuZfi = c(« ) ^ 



q 

This quantity reduces to the linear combination of (qGq) condensate contributions to the correlator for vector currents 
(see (0,|))- In this way, we get the formula 



2 



-320(g-z) 
U (qGq) = jr-g na s {qq 



\2 



81q 6 

from which we then obtain the fraction 80/81 appearing in Eqs. (14)-(|l5|). 

If SD Ml/;Q; g(g) = 0, then we would have the same SRs both for p- and &i-mesons. We process this hypothetical SR 
within the standard approach without a s -corrections in the perturbative contribution and obtain the following values 
for the low-energy parameters of a hypothetical p&i-meson in Anti-Dual Nature: 



(0.865 ± 0.030) GeV ; f pbl = (0.162 ± 0.005) GeV ; s pbl = 1.58 GeV 2 . (31) 



We see that the mass and the decay constant of the p-meson are not so much affected by this (anti-)duality breakdown 
(10% for the mass). The case of the 6i-meson is quite opposite. Here the mass falls down to 45% (the decay constant, 
to 16%, see (f2~ij)). This seems to be quite natural. In the case of the p-meson, the deformation of the SR is large (the 
quark condensate contribution is enhanced by a factor of 3.25), but its functional dependence on the Borel parameter 
M 2 is almost the same. This is not the case for the &i-meson. The deformation of the SR due to the opposite sign of 
the quark condensate contribution is essential and this results in such a large effect for the mass of the 6i-meson. 



QCD VACUUM TENSOR SUSCEPTIBILITY 



The QCD vacuum tensor susceptibility \ has been introduced in Js],^) in order to analyze, in the QCD SR approach, 
the nucleon tensor charges g% and g^- It is defined through the correlator ([j]) 
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X=^, U X (q 2 )^U^(q 2 ) . (32) 
He and Ji || obtained for n x (0) the value 

Y^ n x(°) ~ °- 002 Gey2 ■ ( 33 ) 

In order to obtain a reliable estimate in our approach we substitute the decomposition (|^) in (|32"|), use the relation 
( |30| ) and arrive at the expression 

n x ( g 2 ) = 3(n + ( g 2 ) + n_( (? 2 ))-6n SD ( (? 2 ) . (34) 

This relation clearly demonstrates that IT X is formed by the SD part of OPE, i.e. , by the four-quark condensate 
contribution. Using the dispersion relations (^) 

i2 nx(0) = I s ds ' (35) 

and the phenomenological models for spectral densities p±(s) in (j|), the value of n x (0) can be expressed in terms of 
mesonic static characteristics (the analogous formula has been published in 0]), 



1 Hx(0) = (fD 2 ~ (/ P T ) 2 - (J?) 2 + W - s bl = I -0.0055 ± 0.0008 GeV 2 [NLC] ^ 



'pi \Jjsj_ , J p,p' 

12 2 ,(>;t2 1 -0.0053 ±0.0021 GeV 2 [Standard] 



presented in (p^|)~jl9|) for NLC SR and (pfj|)-(|2l]) for the standard SR respectively^. Note here, both the resuts are 
very close one to another due to strong cancellations in the difference (|3^). So, just this combination accumulates the 
effect of the four-quark part of the whole condensate contribution. If we return to the example of an anti-dual Nature 
(see the end of the previous section, (|3l|)) where this contribution is absent, we obtain the exact cancellation in fl3C|), 
i.e. II x (0) = 0. 

B&O in (|] have used the specific representation that leads to the decomposition 

n x (g 2 ) = 12n 1 ( 9 2 ) + 6g 2 n 2 (g 2 ) , (37) 

where n^g 2 ) = ill + (g 2 ) and q 2 U 2 (q 2 ) = \ (n_(q 2 ) -Tl + (q 2 )). Erroneously suggesting that lim [<7 2 n 2 (<7 2 )l = 
and using the trick suggested in |l^] based on the dispersion relation^] 

1 1 r°° o phcn (s) - o pert (s) 

n x -p-(o) ^ in x (o) = i / £ — ( ' s) p (s) d s , (38) 

they concluded that 

H x (0) n -P- = n?' p '(0) = {flf - J^f « -0-008 GeV 2 . (39) 



5 Depicted errors are obtained by a special invented x 2 -criterium and take into account only the SR stabiliy. 

6 Here pP^ en (s) and p pert (s) are the corresponding spectral densities; the difference of these functions validates the usage of 

unsubtracted dispersion relation 
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But we see from our analysis that the value of (n_(<7 2 ) — H + (q 2 )^ is identically equal to only in an absolutely 
self-dual world, which is definitely not realized in QCD 

lM ° ) 2 n+(0) = Sp ' P '^ 2 Sbl (/J) 2 - (#) 2 - Ulf = -0.0060 ± 0.0017 GeV 2 [NLC] . (40) 

This value is comparable with the value of the B&O estimate ( p9| ) for ll"' p '(0) and should be definitely taken into 
account. Comparing two estimates, our ( |36"| ) and B&O (|3^), one sees not so large deviation from one another. It should 
not be taken by surprise because radiative corrections significantly reduce the B&O value to n°' p '(0) rs —0.003 GeV 2 . 
For this reason the actual magnitude of our total correction to this estimate is of an order of 100%. When our paper 
was finished we find the paper 01 which contains an estimate of the correlator, ^n x (0) = -(0.0083-0.0104) GeV 2 , 
using the constituent quark model. The authors of the paper have determined also a rather wide window for VTS 
by analog of Eq.@ using QCD SR results from different literature sources: jgll x (0) = -(0.0042 - 0.0104) GeV 2 . 
As we pointed out in the Introduction, since these meson constants appear in VTS in a form of a difference, one has 
to define them more precisely and in the framework of a unified approach. So the large width of this window is not 
surprise for us. 

Finally, let us briefly discuss the effect of our estimate of VTS on the nucleon tensor charges. Here we follow to 
pioneering paper by He and Ji Q where these charges were roughly estimated using two types of SRs. Our result ( [36] ) 
increases the lower (decreases the upper) boundary for the g^ {gj) charge approximately by a factor of 1.4: 

$ = 1.47 ±0.76; (41) 
g% = 0.025 ± 0.008 . (42) 

(The results of He and Ji gj, = 1.33 ± 0.53 and g^ — 0.04 ± 0.02 have been obtained for too low, in our opinion, value 
of A QCD = lOOMeV. We, instead, use the value of A QCD = 250MeV.) 
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Appendix 

APPENDIX A: EXPRESSIONS FOR NONLOCAL CONTRIBUTIONS TO SR 

The form of contributions of NLCs to OPE ([ll]) depends on a model of NLC. At the same time the final results 
of SR processing demonstrate stability to the variations of the NLC model provided the scale of the average vacuum 
quark virtuality A 2 is fixed. Here we use the model (delta-ansatz) suggested in jjj and used extensively in [gj; the 
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model leads to Gaussian decay for scalar quark condencate [], (q(0)E(0, z)q(z)) ~ (qq) exp (— |,z 2 |A 2 /8) (see details 
in dominated in b± via 64. Here the factorization hypotheses is applied for the four-quark condensate reducing its 
contribution to a pair of scalar condensates. In the NLC approach this may leads to an overestimate of the four-quark 
condensate contribution due to an evident neglecting of a correlation between these scalar condensates, see Fig. 2. 
In this model we obtain the "coefficients" for OPE in the SR @-(|l7|) 

h^p. = h (M 2 ) + b 3 (Af 2 ) ± 6 4 (Af 2 ) , (A.l) 

where 6 2 corresponds to the vector {(q^mq)), 63 to the quark-gluon-quark ((qG^ v q)) and fe 4 to the four-quark ((qqqq)) 
vacuum condensate contributions (here A = A 2 /(2M 2 )), 

b 2 (M 2 ) = -16 ; (A.2) 
, 4 (60 - 273A + 359A 2 - 134A 3 ) , 

i3(M)=J mT^? [ (A ' 3) 



log (1 - A) ^5A 2 - 21A - 6 
A + (1 - A) 2 



64 (M 2 ) = 24 (A - 7) + 4 — . (A.4) 



The gluonic contribution a± coincides in this model with the standard expression (|13|). For quark and gluon conden- 



sates we use the standard estimates (for "renorm-invariant" quantities in (A. 5) we do not refer to any normalization 
point) 
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! -G 2 ) = 1.2 ■ 10~ 2 GeV 4 , (v^Igg) 2 = 1.83 • 10~ 4 GeV 6 , (A.5) 



A 2 U 2 « 1 GeV 2 ) = (jftggg^g) = 0.4 ± 0.1 GeV 2 . (A.6) 

1^ 1 (qq) 2(qq) V ' 



APPENDIX B: INPUT PARAMETERS IN B&B AND B&O PAPERS 

Both groups of the authors of || and Q used different definitions of initial parameters for processing the SRs. 
Namely, B&O used the following set of values (without any indication on the scale at which renormalization non- 
invariant quantities are determined): 

a s w O.Itt = 0.314 , A-K 2 (qq) = —0.55 GeV 3 , (^/a~ s qq) 2 w 0.61 • 10~ 4 GeV 6 , (B.l) 

whereas B&B0 (on the scale [i 2 « 1 GeV 2 ) 

a s = 0.56, (qq) = (-0.250) 3 GeV 3 , (y/c^qq) 2 ~ 1.37 • 10~ 4 GeV 6 . (B.2) 

This resulted in approximately the same overall contributions of the quark condensate in both papers, see the end of 
Sect. 2. 



7 Here E(0,z) — Pexp(i J dt^A^^Ta) is the Schwinger phase factor required for gauge invariance. 
8 Let us remind that the standard value is ( v / o7gg) 2 ~ 1.83 • 10~ 4 GeV 6 Q. 
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